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The self-energy, spectral functions and susceptibilities of 
2D systems with strong ferromagnetic fluctuations are consid- 
ered within the quasistatic approach. The self-energy at low 
temperatures T has a non-Fermi liquid form in the energy win- 
dow < Ao near the Fermi level, where Ao is the ground- 
state spin splitting for magnetically ordered ground state, and 
Ao oc T^^'^ \n^^'^ {vf /T) in the quantum critical regime («f 
is the Fermi velocity). Spectral functions have a two-peak 
structure at finite T above the magnetically ordered ground 
state, which implies quasi-splitting of the Fermi surface in the 
paramagnetic phase in the presence of strong ferromagnetic 
fluctuations. The triplet pairing amplitude in the quasistatic 
approximation increases with increasing correlation length; at 
low temperatures T <^ Aq the vertex corrections become im- 
portant and the Eliashberg approach is not justified. The 
results for the spectral properties and susceptibilities in the 
quantum critical regime near charge- (spin-) instabilities with 
large enough correlation length ^ ^ {T/vf)~^^^ are obtained. 



I. INTRODUCTION 

Anomalous non-Fcrmi-liquid behavior of correlated 
low-dimensional electron systems has attracted much at- 
tention during the last decade. This behavior is usu- 
ally connected with the violation of the quasiparticle (qp) 
concept in some energy window around the Fermi level. 
A prominent example is the pseudogap phenomenon ob- 
served in underdoped high-Tc compounds [1]. While 
antiferromagnetic (AFM) fluctuations may be responsi- 
ble for non-Fermi-liquid behavior and superconducting 
pairing in cuprates, there is number of systems where 
ferromagnetic (FM) fluctuations may play an important 
role. In particular, FM fluctuations may be important in 
some triplet superconductors, such as UGe2, ZrZn2 and 
Sr2Ru04. These systems motivate studies of electronic 
properties in the vicinity of a FM instability and their 
influence on the triplet superconductivity. 

Although many results exist for the electronic prop- 
erties in the vicinity of an AFM state [2-8], much less 
is known about the evolution of quasiparticle proper- 
ties near the FM instability. An e^/'^ energy dependence 
of the self-energy at the quantum critical point (QCP) 
can be derived from calculations in the context of gauge 
field theories [9], the phase separation problem [10], and 
the Pomeranchuk instability [11], which are expected to 
have the same structure of self-energy corrections as a 
FM instability. The breakdown of the qp concept at the 
QCP is even more apparent at finite temperatures. It 



was demonstrated for fermions interacting with a gauge 
field that the imaginary part of the self-energy in a non- 
self-consistent calculation diverges at the Fermi level at 
T > as a consequence of the divergence of the gauge 
field propagator at zero momentum and frequency [9]. 
Similar behavior induced by the divergence of the static 
uniform spin susceptibility x(0, 0) can be expected for the 
zero-momentum particle-hole instabilities of fermion sys- 
tems with short-range interactions. This behavior can be 
especially pronounced in the renormalized classical (RC) 
regime [12], where the correlation length ^ is exponen- 
tially large. 

The self-energy and the spectral functions in the RC 
regime in the vicinity of a FM instability were previously 
studied within the two-particle self-consistent (TPSC) 
approximation [13,14], one-loop functional renormaliza- 
tion group (fRG), and Ward-identity approaches [14]. It 
was argued, that spectral functions have two-peak struc- 
ture analogously to the vicinity of an AFM instability 
[6] . Contrary to the situation in the vicinity of an AFM 
instability, however, the abovementioned two-peak struc- 
ture of the spectral functions does not imply strong sup- 
pression of the density of states at the Fermi level, but 
leads to the quasi-splitting of the Fermi surface at low T 
already in the paramagnetic phase [14]. While the treat- 
ment based on the TPSC and one-loop fRG approaches 
does not account for the feedback of the self-energy ef- 
fects, the analysis of the self-energy and vertex correc- 
tions using Ward identities has shown [14] that these two 
types of effects almost cancel each other, and therefore 
resulting spectral functions closely resemble their form in 
non-self-consistent approaches. 

These anomalous spectral properties may have a pro- 
found effect on the triplet superconductivity. One can 
expect that due to strong FM fluctuations in the RC 
regime the triplet pairing will be mostly enhanced at the 
new preformed Fermi surfaces. Anomalous spectral prop- 
erties may have important influence also in quantum- 
critical (QC) regime [12], where the quasi-splitting of 
the FS is absent. Previous investigations of this regime 
[15-17] neglected vertex corrections, which may be im- 
portant for large enough correlation length. 

To consider anomalies of electronic properties and their 
impact on the triplet superconductivity, one needs a tool 
which is able to consider self-energy and vertex cor- 
rections on the same foot. The abovementioned self- 
consistent treatment of self-energy and vertex corrections 
near a FM instability was performed only to flrst order in 
1/M, M being the number of spin components (Af ~ 3 
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for the Hubbard model) . It appears important to investi- 
gate spectral properties in the vicinity of FM instability 
beyond the leading order in 1/M to verify whether the 
near-cancellation of the self-energy and vertex corrections 
persists also in higher orders of the 1/M expansion and 
to investigate the effect of the anomalous properties on 
the triplet superconductivity. 

Due to an almost static character of spin fluctuations 
at large correlation length, a useful nonperturbative tool 
for calculation of the spectral properties and susceptibil- 
ities in this case is the quasistatic approach. This ap- 
proach was originally proposed for the summation of di- 
agrammatic series for the self-energy of one-dimensional 
(ID) systems in the vicinity of the charge-density wave 
instability [18] and further developed for 2D systems in 
the vicinity of an AFM instability [3,19]. The quasistatic 
approach allows to sum up the most important series 
of static contributions to the self-energy and interac- 
tion vertices. This approach becomes exact in the limit 
^ ^ oo, and can be applied to study spectral properties 
in the RC and QC regime, provided that the correla- 
tion length is sufficiently large, {(,/a)~^ <C {T/vfY^^, 
where vp is the Fermi velocity, a is the lattice constant 
(the latter criterion follows from the comparison of static 
contributions to the scattering rate ~ TS^/a with the 
dynamic contributions proportional to Vf{T /vpY^^ , cf. 
Ref. [14]). Although the quasistatic approach was applied 
previously to systems in the vicinity of a FM instability 
[20], only the form of spectral functions was analyzed, the 
self-energy, magnetic and triplet pairing susceptibilities 
being not investigated. 

In the present paper we apply the quasistatic approach 
to 2D systems with nonsingular density of states, which 
are on the verge of a ferromagnetic instability, to study 
spectral properties and the possibility of triplet pairing 
in these systems. In Sect. II wc concentrate on the ana- 
lytical results for spectral properties and susceptibilities 
for linear electronic dispersion at ^ — > oo and compare 
these results with the results at finite correlation length. 
In Sect. Ill wc consider the two-particle properties: mag- 
netic susceptibility and the susceptibility with respect to 
triplet pairing. In Sect. IV we discuss the application 
of the results to the quantum-critical regime. Finally, in 
Sect. V we summarize main results of the paper. 



II. SPECTRAL PROPERTIES IN THE VICINITY 
OF FERROMAGNETIC INSTABILITY 



We consider a spin-fermion model [2,3] 
Z[in\= J D[c,c^]D[S]exp{-S[c,S,Tj]) 

S[C, S,n]=J2 [cLii^n - £k)Cka - {ci^Vka + vla^ka) 

k 

+T-i^(x-i + c/^n,)(s,s_,) 



+ UT-^y2^k-k'(Taa'cl^Ck'o 



kk 



(1) 



where q = (q, «w„) and similar for k, LOn = (2n + 1)7tT 
are fermionic Matsubara frequencies, £k is the electronic 

dispersion, a^a-i arc Pauli matrices, Xg — nq/(l — UHq) 
is the dynamical spin susceptibility. 



/(ek+q) - /(£k) 



£k+c 



■ek 
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is the bare polarization operator, /(e) is the Fermi dis- 
tribution function, U is the strength of the interaction 
of electrons with the collective magnetic excitations, the 
lattice constant a = 1. Although this model was orig- 
inally proposed as phenomenological model for systems 
with strong AFM fluctuations, it can be applied for sys- 
tems with strong FM fluctuations as well. Generally 
speaking, the interaction U differs from the bare on-site 
Coulomb repulsion because of contributions of the chan- 
nels of electron-electron scattering different from particle- 
hole one, therefore U should be considered as an effec- 
tive interaction. The counterterm proportional to U^Ilq 
keeps the renormalized spin-spin propagator equal to Xg 
(see below). The rigorous derivation of the model (1) from 
the microscopic Hubbard model will be considered else- 
where [22]. 

Integrating out fermions from (1), we obtain 

Z[r^]^ J D[S]eM-Seff[S,v]) (3) 
Seff[S,v] = T-i ^x,-'S,S_, - lndet[-G-,l(S)] 



kk' 



where 



^kk'i^) — (^'^" ~ Sk.)S<T<T'^kk' + USk-k'O-o 



(4) 



In the following wc expand lndet(— G^^) in Eq. (3) in 
powers of S and retain only quadratic term, which is 
exactly cancelled by the counterterm introduced in Eq. 
(1) (so that the propagator of the field S remains equal 
to Xg), the relevance of higher-order terms is discussed 
below. 

We start investigation of the functional (3) with the 
consideration of the limit ^ — > oo, where spin fluctua- 
tions are especially strong. At ^ ^ oo the susceptibility 
Xq = x(q;iu;„) is divergent at q = 0, ujn = 0. Since the 
momentum transfer for the scattering on these most sin- 
gular magnetic fluctuations is small, it can be neglected 
in the electronic Gcen functions. Sums over internal mo- 
menta in all diagrams are applied then only to the prop- 
agators Xq of the spin field Sq, so that the action (3) in 
^ ^ oo limit can be reduced to an effective action which 
contains only one fluctuating field S = 85=0 (cf- Refs. 
[18,3,19]) 
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Z[n]= J d^SeM-Seff[S,v]) 



(5) 



t f iujn-ek + US'' 



US+ 



us- 



where = ± iS^. The effective propagator of the 
field S, Aq/{3U^), is determined by the average (local) 
spin susceptibility 



(6) 



For an ordered ground state Aq is almost independent 
of temperature at low T and its T ^ limit is equal 
to the square of the ground-state spin splitting. At the 
same time, in the QC regime we have Aq oc T\n{vF/T), 
in this case the effect of finite correlation length should 
be also accounted for. Due to neglection of terms which 
are of higher order n > 2 in spin operators and propor- 
tional to Pq" ^'(c), [/5o(e) being noninteracting density 
of states], the generating functional (5) is valid only for 
regular Po{e), which is smooth enough in the vicinity of 

the Fermi level to satisfy U"'pQ^\e) <C po{£) at n > 
[22]. 

Similar to Ref. [3] we generalize the action (5) to M- 

component field S = {Si...Sm), M = 3 for the model (1); 
this generalization also allows one to consider a charge 
instability with M ^ 1. The results for the observable 
quantities are found by differentiation of partition func- 
tion over the source fields rj and are expressed as integrals 
over the field S of some functions f{S). 

For the electronic Green function at ^ — > oo we obtain 
the result 



G{cJ) = 



= j d'^S 



1 / Maf 

W \ 2A2 



exp(-5e//[S,0]) 



M/2 



M 

-or(i-i^. 



(7) 



which depends on cU = w — ffk only, F (a, x) is the in- 
complete Gamma function. The result (7) is similar to 
previous result in the vicinity of an AFM instability [3]. 
The electronic self-energy is given by 



(8) 



The retarted Green function and self-energy on the real 
axis are obtained by the replacement oj ^ u> + iO~^. For 
the following analysis it is convenient to introduce a one- 
particle irreducible (IPI) vertex 



7''(k,zu;) = G^ 



k'(T"cr" 



/3 rf3 

dr / dr'e' 
^0 



(9) 



X limx-^(r[4,,.(T)ck'+q,.'"(T')S«(0)]>iPi 



where a = x,y,z. Similar to the Green function (7), 
7(k, w) depends on ZJ only: 7"(k, cj) = jiZJ)- The func- 
tion 7(tJ) can be obtained from the exact Dyson relation, 
connecting the vertex and the self-energy (see, e.g., Ref. 
[23]), which at ^ ^ 00 takes rather simple form 



S(a;) 



Ag7(aJ) 
u) - E(ZZJ) 



(10) 



The quantities G{u)), S(a;), and 7(0;) determined by Eqs. 
(7)-(10) can be considered as perturbation series in Aq oc 
U. The corresponding lowest-order coefficients obtained 
from Eqs. (7) and (8) are 



S(a;) 
7(0;) 



^0 J 

CO Mw3 M2w5 ^ 

(2 - M) A§ 2 (4 - M) A^ 



2A^ , 2(4 + M)A8^^^^,^ 



1 + 



+ 



+ 0{Al) (11) 



The coefficients of the series in Aq can be found also 
directly from a diagram technique (we have verified cor- 
respondence of several lowest-order terms) . 

The perturbation scries (11) breaks down at frequen- 
cies |a;| < Aq, although nonperturbative results (7)-(10) 
can be used to analyze physical properties in this fre- 
quency range as well. In particular, for M > 2 we find 



ReS(w) ~ (M - 2)Al/Muj, \uj\ « Aq. 



(12) 



Comparing this result with the perturbation theory result 
(11) one observes, that ReS(a;) has a nontrivial crossover 
with the reduction of the number of spin components at 
low frequencies. Such a crossover is similar to that for the 
spin-spin correlation function in 2D and quasi-2D gener- 
alized Heisenberg model with 0{M)/0{M — 1) symmetry 
[24,25]. In the crossover region |a;| ~ Aq the real part of 
the self-energy S(ti') is only weakly ij-dependent. The 
imaginary part of the self-energy at small |a;| and M > 2 
reads 



ImE(w) ~ -Am(2/M)1/2(m/2 - l)^Ao 
x(Ma;V2A2)(^-3)/2 
-7r(M - 2)A^^(a;)/M, jwj -C Aq 



(13) 



where Am = n/{M/2 - 1)1 for even M and Am = T{1 - 
M/2) for odd M. For the charge instability case (M = 1) 

and small |w| <C Ao we obtain 



ReS(a;) ~ 2a;/7r2; ImS(a;) 



{M = l) 



The overall frequency dependence of £(0)), 7(0;), and the 
spectral function A{cl)) = — ImG(a;)/7r at ^ — > cxd calcu- 
lated using Eqs. (7)-(10) for different M is shown and 
compared with the results of 1/M expansion of Ref. [14] 
in Fig. 1. One can see that at M = 3 (the same behavior 
takes place for all M > 1) the real part of the self-energy 
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FIG. 1. The real and imaginary parts of the self-energy 
(a,b), the spectral function (c), and the vertex function 7 
(d) in the quasistatic approximation at ^ — > cxd, M = 1 
(dot-dashed lines) and M — 3 (solid lines) as a function of 
ZJ = uj — Ek for Aq/vf = 0.1, compared to the result of 
1/M expansion (dotted lines) and the one- loop functional 
RG approach for spin-fermion model (dot-dot dashed line) 
for M = 3. 



has (infinite) positive slope at the Fermi level, where the 
imaginary part of self-energy has (5-like singularity, and 
the spectral function has two-peak structure. These fea- 
tures are in qualitative agreement with previous results of 
the first-order 1/M analysis [14]. They arise as a result of 
strong FM fluctuations and violate the quasiparticle con- 
cept near the Fermi level. It was argued in Ref. [14] that 
the two-peak structure of the spectral function, together 
with its dependence on w — £k implies pre-formation of 
the two new Fermi surfaces already in the paramagnetic 
phase (so called quasi-splitting of the FS), the same ar- 
guments can be applied to the results of quasistatic ap- 
proach. 

The main difference of the results of quasistatic ap- 
proach from the results of 1/M expansion [14] is in par- 
tial transfer of the spectral weight from the peaks of the 
two peak structure to small w-region, where the spectral 
weight in the result (7) is small, but finite. From Eq. (7) 
we find A{uj) ^ at small w. The nonanalytical 

dependence of A{uj) on M explains why this behavior is 
not captured by 1/M expansion. 

For M = 1 (charge instability case) the imaginary part 
of the self-energy is finite at the Fermi level [see Eq. 
(13)] and the spectral function has one-peak structure. 
It does not imply, however, the validity of the quasipar- 
ticle concept, since the real part of the self-energy has 
positive slope at the Fermi level, which invalidates quasi- 
particle picture. Note that vertex corrections are finite 
in this case, and, therefore, are not as important, as for 



M > 1. Since the long-range order exists for M = 1 also 
at finite T, these results are applicable only in a narrow 
critical regime near the transition temperature, but have 
also some implication for quantum-critical region, as dis- 
cussed in Sect. IV. The behavior of G{uj), S(w) and 7(w) 
at M ~ 2 (XY-type symmetry) is very similar to that for 
M = 3, except for additional logarithmic corrections. 

It is instructive to compare the results (7)-(10) for the 
self-energy and vertex at ^ ^ cxd with the corresponding 
results of recently proposed functional renormalization- 
group approach for the boson- fermion model [26]. Since 
the momenta integrations and frequency summations in 
Feynman diagrams are restricted at large ^ to iiOn — 
and the near vicinity of g = 0, neither frequency, nor 
momentum cutoff of electronic or bosonic degrees of free- 
dom are convenient for this problem. Instead, we impose 
the temperature cutoff on the electronic Green function 
(the correlation length is kept fixed, so that the bosonic 
propagator does not acquire temperature dependence). 
We also combine this scheme with the one-particle self- 
consistent modification of fRG equations [27], which al- 
lows for a correct treatment of self-energy effects to one- 
loop order. The resulting one-loop fRG equations at 
^ — ^ 00 read 



duj 
duj 



M 



did 



(14) 



dG{uj) 
doj 



where G{uj) = [w — E(w)]^^. We compare the solution of 
Eqs. (14) with the result of quasistatic approach (7) in 
Fig. 1 . One can see that the one- loop fRG equations de- 
scribe very accurately the perturbative regime |a;| > Ag, 
but their description breaks down in the strong-coupling 
regime |a;| < Aq. 

To study the effect of finite correlation length, we em- 
ploy an ansatz for the nonuniform magnetic susceptibility 



x(q,o) = 



A 



(15) 



Note that this ansatz neglects recently found nonanalytic 
corrections [21] and is, therefore, valid above the charac- 
teristic temperature Tx ~ U'^/vp, where these correc- 
tions become important. At finite ^ the quasistatic ap- 
proach can be applied when static contributions to the 
self-energy and vertices are dominating near the Fermi 
level, i.e. at ^ {T/vpY^^, as discussed in the intro- 
duction. This condition is satisfied, in particular, in the 
RC regime. 

The generalization of the quasistatic approach to the 
susceptibility ansatz (15) is considered in Appendix. The 
result (7) is to be replaced at finite ^ by an integral recur- 
sion relations for the electronic self-energy S(u;) = Si(w) 
and vertex 7(0;) = 71 (w). 
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FIG. 2. The real and imaginary parts of the self-energy 
(a,b), the spectral function (c), and the vertex function 7 (d) 
in the quasistatic approximation at M = 3, Aq/vf = 0.1, and 
different values of the correlation length. 



7j(w) = 1 



Agr, da 



(17) 



where Gj{uj) = [u — 



Cj = j'/M (even j), 



c^- = (j 4- M - l)/M (odd j); = j/(M - 2) (even 
j), r, = (M ~ 2)(j + M- l)/AP (odd j), Goo(cc') = l/o;, 
and 7oo(w) = 1. The most important contributions to 
integrals in Eqs. (16) and (17) at ^ — > 00 come from a 
narrow vicinity of the point a = 0, and these equations 
reduce to the continuous fraction representation of the 
gamma- function in Eq. (7). At the same time, at finite 
f the Eqs. (16) and (17) have to be solved numericaUy. 

The results for the self-energy and the vertex are shown 
for different values of ^ and compared with the results 
for ^ 00 in Fig. 2. In agreement with previous anal- 
ysis [14], the real part of the self-energy acquires a large 
positive slope at the Fermi level, dKeYi/de ~ T^^. The 
(5-function singularity at ^ ^ 00 in the imaginary part 
of the self-energy is replaced by the lorentian-like form 
of the imaginary part with |lml](0)| ~ T^, so that the 
quasiparticle picture is invalid at finite ^ as well. With 
decreasing correlation length the structure of the spec- 
tral function changes from the two- to one-peak form at 
~ Aq/vf- Contrary to ^ — > 00 limit, the vertex re- 
mains finite at finite ^. 

At M = 1 the imaginary part of the self-energy, which 
was finite at ^ — > 00, is determined by |ImE(0)| ~ 
min(Ao,T^). The qp picture is violated in this case as 
well, since the slope of ReE is positive, dRe'S/de ^ 1. 



III. TWO-PARTICLE PROPERTIES 

Now we discuss two-particle properties. First we con- 
sider static uniform spin susceptibility (/>ph. According 
to Ref. [23], this susceptibility can be expressed through 
the irreducible in the particle-hole channel susceptibility 
0ph,o via the relation 



0ph = 't'phfi/ {i-U (j)ph,o) 



(18) 



which is similar to the random phase approximation re- 
sult with the difference that (pph.a includes the self-energy 
and vertex corrections. Using the definition of the irre- 
ducible vertex, Eq. (9), we find 



0ph,o = J depo{e)^{e) 

$(e)--T^0(zc.„-e), <j>{z) ^ G'{zh{z) (19) 

The necessary condition for the existence of a ferromag- 
netic instability is 0ph,o > 0; the function $(e) character- 
izes the relative weight of states with different energy in 
9^ph,o- The plot of the function 0(z) in the complex plane 
at ^ — > oo is shown in Fig. 3a (the plot of this function 
at finite ^ looks similarly). At ^ ^ cxd the contribution 
of regions < Ao and \uj\ > Aq to ^(e) have differ- 
ent signs and compensate each other for all e, except for 
[el < Aq, where $(e) is maximum (Fig. 3b). Therefore, 
irreducible susceptibility depends on the details of the 
density of states only in the energy range |e| < Aq. One 
can see that for regular densities of states, which are not 
strongly suppressed in this energy range, the condition 
0ph,o > can be easily fulfilled. Stronger criterium of 
stability of ferromagnetism d'^(f)q/dq^ < is studied in 
detail in a forthcoming paper [22]. 

To investigate the static magnetic susceptibility at fi- 
nite ^, we suppose that the temperature dependence of 
the correlation length is given by ^ = exp(r*/r) where 
T* is the crossover temperature to the renormalized clas- 
sical regime. The function $(£) for different values of ^ 
is shown in Fig. 3b. At not too large correlation length 
the energy range which contributes to <^ph,o is spread to 
[el > Aq as well. At the same time, the total area un- 
der <i>(e) changes rather weakly and, therefore, one can 
expect weak dependence of the irreducible susceptibility 
on the correlation length ^. 

The susceptibility with respect to triplet pairing 



^PP.tr = E E K.'A^,,,, f dr 

11/ ////// -'0 
kk' (T(T'a"(T"' 

x(T^[cLWclk,.'Mck',."(0)c_k' 



(20) 
'(0)]) 



where = [a°'ay)aa' can be considered in a similar 

way. It is convenient to represent it in the form 



^pp,tr 



tE E<-"^;pTr'(k,^) (21) 



5 



b) 6 






b) 



800 



FIG. 3. (a) The plot of the function cj){iu! — e) at ^ — > oo 
and Ao = O.lvp- (b) The function ^(e) (which determines 
the uniform static irreducible spin susceptibility according to 
Eq. (19)) at Ao = O.lvp, T* = 0.25i;f, and different values 
of the correlation length. 
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Iz.' .a" a'' 

xc_k',."'(0)]) (22) 

Similar to magnetic susceptibility, 4'p'ptr depends on 
and uj only and can be generally written as 

(/•ppTr (k, it^n) = A;^^,0pp,tr(ek, iuj) (23) 
0pp,tr(ek,«t^) = G{iuj - eu.)G{-iujn - ek)7tr(ek,«w) 

where jtr{sii^,iLu) is the pairing vertex. At ^ ^ oo we 
obtain 



0pp,tr(£k, »w) = F(ek, t^) + F{-e-k, u) 



(24) 



where 



uj+{M -l)sk 
2 M Lj ZUeic 



1 - 



M 



2A2 



(25) 



Note that jtr depends on lu and separately and 
7tr(0,a') = 7(0;). The pairing susceptibility (20) is ex- 
pressed through 0pp,tr(£k, *w) by the relation 



$tr(£) = (/)pp,tr(£, iuJn), 



The function 0pp.tr(£, iw), which characterises the contri- 
bution of different momenta and energies in the pairing 
susceptibility is plotted at ^ — > oo in Fig. 4a. This func- 
tion is divergent at w ^ 0^ , which signals the possibility 
of the triplet pairing at T 0. At finite small uj the 
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FIG. 4. (a) The plot of the function (/!>pp,tr(e, itij) at ^ oo 
and Ao — O.lvF (b) The function 0pp,tr(£, ittT) at finite ^, 
Ao = O.l/S^/^UF, and T* = 0.05vf. (c,d) The function 
"I>tr(e), which determines the triplet pairing susceptibility ac- 
cording to Eq. (26), at the same values of Ao and T* as in 
b) (c) and at Ao = Q.lvp, T* = 0.25uf (d). 



function 0pp_tr(£, iti^) is maximum at £ = ±Ao. There- 
fore, contrary to standard BCS problem, the pairing due 
to the ground-state FM instability in 2D system involves 
particles with finite energy (with respect to the para- 
magnetic Fermi surface) and the momenta at the new 
preformed Fermi surfaces. 

At finite ^ the function 0pp.tr (£k, ii-^) is determined by 
the Eq. (23); the pairing vertex 7tr is obtained from 
the recursion relation which is similar to the recursion 
relation for 7, 



7tr,i(£k,^) = 1 



A 



0' 1 



21ne 



da 



X7tr,i+i(£k - VFa,uj)G.j+i{uj 
xGj+i(— a; — £k — vpa) 



£k 



vpa) 
(27) 



with 7tr,oo(t^,£k) = 1- 

The function (/ipp tr(£, ^ttT) at different values of ^ is 
shown in Fig. 4b. With decreasing ^ the two-peak struc- 
(26) ture of 0pp,tr(£, ittT) continuously changes to a one-peak 
structure at T = T*/ln^ ~ Aq. The function $tr(£) 
for two choices of T* and Aq and different values of the 
correlation length is shown in Figs. 4c, d [we rescale the 
value of Ao oc (T*)^/^, as it follows from Eq. (6)]. Similar 
to (/)pp,tr(£, iT^T), the function $tr(£) changes its behavior 
from the two-peak to a one-peak structure at T ^ Aq, 
so that the triplet pairing fluctuations are dominating at 
not too large ^ at the paramagnetic Fermi surface. 
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FIG. 5. The triplet pairing vertex 7tr(e,j7rr) in the qua- 
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sistatic approach at Aq = 0.1/5 
Ao = O.lvF, T* = Q.2'bVF (b), and different 5. 



To clarify the role of the vertex corrections for the 
triplet pairing, we plot in Fig. 5 the triplet pairing ver- 
tex 7tr(e, ittT) at first Matsubara frequency for the same 
choices of T* and Aq as in Fig. 4. We find that at 
T* <C vf (weakly FM ground state) and sufficiently large 
correlation length In^ » T* /Aq the triplet pairing ver- 
tex 7tr(e, JTrT) is considerably enhanced. In particular, 
we emphasize that the divergence of the pairing suscepti- 
bility at ^ — > oo arises solely from the vertex corrections. 

The triplet pairing vertex 7tr(e, ittT) in the quasistatic 
approach can be furthermore compared with the result 
7t^(e,i7rr) of the approach which accounts for the self- 
energy corrections only (the analogue of the Eliashberg- 
type approach of Refs. [16,17]). At ^ ^ oo we obtain in 
such an approach (cf. Ref. [14]) 

JfAe,^0+) = [l-iM-2)Al\GEie + ^0+)\yM]-' 
Ge{s + iO+) = [e-Y.E{s + iO+) + iO+]-^ 

~ (28) 



S£(e + iO+) = {e- Ve-2AoV£ + 2Ao)/2, 



the corresponding finite-^ result can be obtained from 
Eqs. (16), (27) with Cj = 1 and = {M ~ 2)/M. 
It can be found from Eq. (28) that j^.{e,iO+) < M/2 
and, therefore, remains finite at T ^ 0, ^ ^ cxd. At the 
same time, the triplet pairing vertex in the quasistatic ap- 
proach is divergent in this limit, leading to the divergence 
of the triplet pairing susceptibility. This divergence indi- 
cates the possibility of the triplet pairing due to classical 
spin fluctuations, which is complementary to pairing due 
to quantum spin fluctuations previously studied in Refs. 
[16,17]. At In^ < T*/Aq we find jtrie,inT) ~ 1 and the 
vertex corrections are not important. In this case, the 
Eliashberg-type approach of Refs. [16,17] is justified. 

Therefore, the role of the vertex corrections for the 
triplet pairing depends on the temperature and the value 
of the correlation length, at not too large correlation 
length the vertex corrections can be neglected, but they 
become crucially important at large ^. 



IV. THE QUANTUM-CRITICAL REGIME 

As discussed previously, the static contributions to 
self-energy and vertices dominate over quantum contri- 
butions for sufficiently large correlation length, ^ 
(T/vp)^^^- Provided that this inequality is satisfied, 
one can apply the above consideration to the quantum- 
critical regime as well. As was mentioned in Sect. II, 
in this regime Aq oc T^^^ \n'^'^{vp /T) becomes temper- 
ature dependent itself. The spectral properties in this 
case depend on the value of the exponent a, which de- 
termines the temperature dependence of the correlation 
length according to <^ ^ (T/wf)""- In this respect, 
two regimes can be distinguished: (i) ^ Aq/vf, 
i.e. a > 1/2 and (ii) Aq/vf < < {T/vf)^^^, i.e. 
1/3 < a < 1/2. In the regime (i) spectral functions have 
the two-peak structure, so that similar to the RC regime, 
studied in Sect. II, the Fermi-surface at finite tempera- 
tures is pre-split, while in the regime (ii) spectral func- 
tions have one-peak structure. In both regimes the real 
part of the self-energy has positive slope at the Fermi 
level (9ReI]/9e ^ (T/vf)£,'^, and the imaginary part at 
the Fermi level (associated with the inverse qp lifetime) is 
anomalously large, |ImS| ~ T^, so that the qp picture is 
invalid. According to Millis theory [28], the temperature 
dependence of the correlation length in the QC regime is 
given by oc T^/"^ ln'^^{vF/T), and therefore this type 
of dependence belongs to the regime (ii) . 

For the charge instability case (M = 1) the derivative 
of the real part of self-energy at the Fermi level is positive, 
but finite, 9ReE/9e ^ 1, and the imaginary part at the 
Fermi level |ImE| - Aq ~ T^^^ln^^{vF/T) in the regime 
(i) and |ImE| - in the regime (ii). Note that |ImE| 
in the regime (i) does not depend on the value of the 
exponent a in this case. 

The triplet pairing susceptibility in the QC regime is 
determined by the Eq. (26). The function $tr(£) has a 
one-peak structure similar to that in the RC regime at 
not too large correlation length. According to the results 
of Sect. Ill, vertex corrections to the triplet pairing sus- 
ceptibility are small at T > Aq where the Eliashberg-type 
approach of Refs. [16,17] is justified. The corresponding 
condition in the QC regime coincides, up to logarithmic 
corrections, with the condition of the applicability of the 
susceptibility ansatz (15), T > Tx- At the same time, at 
T <C Tx triplet pairing susceptibility is substantially en- 
hanced over its bare value and vertex corrections can not 
be neglected. The analysis of this case requires consider- 
ation of the nonanalytic corrections to magnetic suscep- 
tibility, which is beyond the scope of present paper. One 
can expect, however, that in this regime magnetic and 
superconducting fluctuations are strongly coupled and 
should be considered on the same foot. 
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V. CONCLUSION 

In the present paper we have studied spectral proper- 
tics and the triplet pairing fluctuations in the vicinity of 
a FM instability. The strong FM fluctuations violate the 
qp concept near the Fermi level, leading to anomalously 
large scattering rate at the Fermi level, |ImS| ^ T^, and 
the positive slope of the real part of the self-energy at 
the Fermi level, dReE/de ~ (T/vp)^'^- Although these 
results coincide with the results of the second-order per- 
turbation theory with respect to coupling of electrons 
with magnetic excitations [14], they take into accoimt 
the self-energy and vertex corrections. Therefore, these 
two type of corrections almost compensate each other, as 
it was concluded before on the basis of 1/M expansion 
[14]. At large enough correlation length :§> Aq/vf 
(Aq is the ground state spin splitting in the RC regime 
and Aq cx T^^"^ in the QC regime) the abovementioncd 
features of the self-energy lead to the two-peak structure 
of the spectral function, which implies the quasisplitting 
of the Fermi surface, as proposed in Ref. [14] . The struc- 
ture of the spectral function changes to a one-peak form 
at < Aq/vf- For M = 1 (charge instability case) the 
spectral function has a one-peak structure at arbitrary 
^. This does not restore the qp picture, however, since 
the slope of ReS remains positive, dUc^/de ~ 1 and the 
imaginary part |ImS| ~ min(Ao, z;^^"^) is finite at low 
T. 

The triplet pairing susceptibility near the FM insta- 
bility is considerably enhanced at low temperatures as 
compared to its bare value. In the RC regime at large 
correlation length the triplet pairing is most strong at 
the newly preformed (quasi-split) Fermi-surfaces, and 
with increasing temperature (i.e. decreasing correlation 
length) the triplet pairing arises at the paramagnetic 
Fermi surface. The vertex corrections to the triplet pair- 
ing susceptibility are important at low enough tempera- 
tures T < inax{Ao,VF^~^}- 

In the quantum critical regime, dynamic contributions 
with nonzero bosonic Matsubara frequency to self-energy 
and vertices can be neglected at low enough temperatures 
for a > 1/3, where the exponent a describes the tempera- 
ture dependence of the correlation length, ^ ~ {T/vf)~"- 
Depending on the value of a, one- or two-peak structure 
of the spectral functions is possible, the former arising for 
1/3 < a < 1/2, the latter for a > 1/2. The qp picture 
is violated for any a < 1. Since, however, the contribu- 
tions of nonzero bosonic Matsubara frequencies are dif- 
ferent only by power of temperature, their contribution 
is expected to be important for a correct quantitative de- 
scription of quantum critical regime. The consideration 
of the triplet pairing shows that the vertex corrections 
in the quantum critical regime can be neglected at not 
too low temperatures and become non-negligible in the 
same temperature range T < U'^/vf, where nonanalytic 
corrections to magnetic susceptibility become important. 
The consideration of this region requires an analysis of 



magnetic and superconducting fluctuations on the same 
foot, which is the subject of future investigations. 

In summary, the quasistatic approximation discussed 
in the present paper allows for a treatment of the self- 
energy and vertex corrections which arise from static 
magnetic fluctuations. In this respect, such an approxi- 
mation has some advantages over 1/M expansion, since 
it does not require M to be sufficiently large. However, 
it can be hardly generalized to include dynamic magnetic 
contributions with nonzero bosonic Matsubara frequen- 
cies. Therefore, a generalization of the 1/M expansion, 
which includes these dynamic contributions, is desirable. 
On the other hand, the generalization of the quasistatic 
approach which includes the effect of van Hove singular- 
ities in the electronic spectrum could provide a possibil- 
ity to describe qualitatively the properties of real low- 
dimensional materials. 



VI. APPENDIX. THE DERIVATION OF THE 
RECURSION RELATIONS AT FINITE 
CORRELATION LENGTH 

In this Appendix we reconsider the extension of the 
quasistatic approach to 2D case when the static magnetic 
susceptibility has the form 



x(q,o) 



(29) 



The early version of quasistatic approach for ID models 
[18] can be directly extended to 2D case only for the 
factorizable form of the susceptibility (cf. Refs. [3,19]) 



(30) 



with and q± being the components of q, parallel and 
perpendicular to the electron momentum k. Although 
the extension of quasistatic approach to the susceptibil- 
ity ansatz (29) was discussed previously in Ref. [3], we 
argue that this extension does not treat correctly loga- 
rithmic corrections, which arise after integration of Eq. 
(29) over q. While these logarithmic corrections are sub- 
leading in the quantum-critical regime, they are crucially 
important in the RC regime, where the correlation length 
is exponentially large. 

To discuss the way of a proper generalization of the 
method, we consider the contribution of a 2A^-th order 
diagram for the self-energy (cf. Ref. [3]) 

S(2^)(k,ic^„)=52iv ^ x(qi,0)...x(qiv,0) 



qi...qN 
2N-1 



(31) 



where g = UT^^^. The coefhcients Rja determine 
whether a-th momentum variable c[a enters j-th elec- 
tronic Green function, see details in Ref. [3]. At large 
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^ ^ 1 the most important contribution to E*^^^^ comes 
from small momenta, and it is sufhcient to expand the 
denominator of Eq. (31) in q. For further convenience, 
we introduce new variables of integration Oa = Qa. cos 9 a, 
where a = 1...A'' {9a is an angle between Qq, and k). The 
integrals over q can be then calculated analytically; using 
the form of the susceptibility (29) we obtain 



for the corresponding noncrossing diagram is different. 
Indeed, substituting the dressed Green function instead 
of the bare one in Eq. (32) with = 1, and taking into 
account that the self-energy depends on w — Ek only, we 
obtain the recursion relation 



S,(a;) 



S(2^)(k^,ia;„) = [Ag^/A^r H 



4tt 



dttN 



f 

J — ( 



da 



^/a?+^-^ 
1 



(36) 



2JV-1 

X n 



N 



a=l 



(32) 



The corresponding result for the factorized susceptibil- 
ity ansatz (30) differs by the replacement + — > 
^(a^ + ^~^) in the denominators of Eq. (32), aa = qa,\\ 
in this case. 

For |a;| vf^~^ one can neglect a„ in the denomina- 
tors of Green functions in Eq. (32) to obtain 

S(2^)(kf,a;) ~ (V/47r)^a;-(2^-i)ln^(Ca;/t;ir), 

> vfC'^ (33) 

To find asympthotic form of the self-energy at small 
\lo\ -C w_f^~^, we shift contours of integrations in Eq. 
(32) to the upper half of the complex plane. The inte- 
grals are then determined by the contributions of branch 
cuts of square roots and 

E(2^')(k^,o.) ^ ^{Ag^/4^re''-'f{{Roc.}), 

H-^vfC^. (34) 

where f{{Rja}) is some function which depends on the 
coefBtients Rja only. One can see, that at \lo\ <C vf£,'~^ 
the self-energy E^^''^^ (k;?, w) does not acquire logarith- 
mic corrections. At the same time, the approach of Ref. 
[3] leads to logarithmic corrections in the self-energy in 
both the limits, ^ vpS,^^ a,nd \lo\ <C vp^^^, due to 
an incorrect factorization of Bessel functions of sums of 
auxiliary variables, used in Ref. [3]. Note, that for the 
ansatz (30), branch cut singularities of the integrands in 
Eq. (32) are replaced by single poles, so that at arbitrary 
\u)\ <^ vf we obtain 



E(2^)(k,a;)^(V/47r)^ [] 



1 



(35) 



with rij 



X^^^^-Rja, which reproduces the result of 
Refs. [3,19]. 

For the form of susceptibility (29) one can develop an 
approximate approach, which becomes exact at |a;| ^ 
Vf^T^ ■ Similar to Refs. [18,19] we approximate the contri- 
bution of any diagram by the contribution of correspond- 
ing noncrossing diagram. Although the multiplicity fac- 
tors are the same, as derived in Ref. [3], the expression 



w — vfo- — Sj+i(a' — vfo) 

where Cj = j /M (even j) and Cj = {j + M — \)/M (odd 
j), Contrary to Ref. [3], this is an integral rather than 
an algebraic relation. The initial condition for Eq. (36) 
is Ecc(t^) = 0, the self-energy is given by S(a;) = Ei(a;). 
For the vertices we obtain similarly 



, . 1 Ag^r, f 



da 



7tr,j(£k, W) = 1 + 



7j+i(a; - vpa) 

[w — Vpa — Ej+i(a; — vpa)]'^ 

Ag^r, da 



(37) 



47r 7-00 Va^+t^ 
X7tr,i+i(ek - vpa,uj)[uj - £k - vpa 
-Ej+i(w - ek - vpa)]~'^[-u) - £k 
—Vpa — 'Sj+i{—u> — su^ — VFa)]~^ (38) 

withrj = j/{M-2) (even j), = {M-2){j+M-l)/M'^ 
(odd j) and 7oo(w) = 7tr,oo(^0 = 1- 

As mentioned above, for the ansatz (30) the replace- 
ment y/al + ^ al^ + £_-^ in Eqs. (36)-(38) should 
be made. The integrals in the Eqs. (36) and (37) can 
be then evaluated analytically, leading to the recursion 
relations of Refs. [3,18]. At the same time, the integrat- 
ing expression of Eq. (38) is nonanalytical in both, upper 
and lower half-plane, and therefore can not be reduced to 
an algebraic form even for the factorizable susceptibility 
ansatz (30). 
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